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Abstract. Darboux transformation is developed to systematically find variable sep- 
aration solutions for the Nizhnik-Novikov-Vcsclov equation. Starting from a seed 
solution with some arbitrary functions, the once Darboux transformation yields the 
variable separable solutions which can be obtained from the truncated Painleve anal- 
ysis and the twice Darboux transformation leads to some new variable separable 
solutions which are the generalization of the known results obtained by using a guess 
i— ^ ■ ansatz to solve the generalized trilinear equation. 
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O ■ To find some exact solutions of integrable systems has focused many mathematicians 

o : 

i . and physicists' attention since the soliton theory came into being at the 1960's. There 

(N ; 

are many important methods to obtain the special solutions of a given soliton equation. 
Some of the most important methods are the inverse scattering transformation (1ST) 
approach]]]], the bilinear form|2[], symmetry reduction ||, Backhand transformation and 

•i-H , 

Darboux transformation []5| etc. In comparison with the linear case, it is known that 
1ST is an extension of the Fourier transformation in the nonlinear case. In addition to 
the Fourier transformation, there is another powerful tool called the variable separation 
method M in the linear case. Recently, some kinds of variable separation approaches 
had been developed to find new exact solutions of nonlinear models, say, the classical 
method, the differential Stackel matrix approach H, the geometrical method M, the 



ansatz-based method ||, RL functional variable separation approach [HI], the derivative 
dependent functional variable separation approach |Tl||, the formal variable separation 
approach (nonlinearization of the Lax pairs or symmetry constraints) [0 and the 



informal variable separation methods |T3|-|r5| 
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Especially, for various (2+l)-dimensional nonlinear physics models, a quite universal 
formula 

jj _ 2(aia 2 - a^z)%Vx ^ 
(ao + dip + d 2 q + d 3 pq) 2 ' 

where a*, (i = 0, 1,2, 3) are arbitrary constants and p = p(x, t) and q = q(y, t) are 
arbitrary constants of the indicated variables, is found by using the informal variable 
separation approach|[L3|-[ 15 . Starting from the universal formula (0]), abundant local- 



ized excitations like the dromions, lumps, ring solitons, breathers, instantons, solitoffs, 
fractal and chaotic patterns are found. Now a very important question is can we find 
the universal formula from other well known powerful methods like the 1ST approach, 
dressing method, Darboux transformation (DT) etc? 

DT is one of the most powerful methods to construct a broad class of considerable 
physical interest and important nonlinear evolution equations such as the well-known 
Korteweg-de Vries (KdV) equation, the Kadomtsev-Petviashvili (KP) equation, the 
Davey-Stewartson (DS) equation, the sine-Gordon (SG) equation (3J and so on. In this 
letter, we use the Darboux transformation to study the variable separable solutions for 



the (2+1) dimensional Nizhnik-Novikov-Veselov system. JO: 



u t = u xxx + u yyy + 3(vu) x + 3(uw) y , (2) 



u x = v y , u y = w x . (3) 

The (2+l)-dimensional NNV equation is an only known isotropic Lax integrable exten- 
sion of the well-known (l+l)-dimensional KdV equation. Many authors have studied 
the solutions of the NNV equation. For example, Boiti et al[|r7| solved the NNV equa- 
tion via the 1ST; Tagami and Hu and Li obtained the soliton-like solutions of the NNV 
equation by means of Backlund transformation [|i~8|j ; Hu also gave out the nonlinear 
superposition formula of the NNV equation[|19j; Some special types of multi-dromion 
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solutions were found by Radha and Lakshmanan[20]; The generalized localized exci- 
tations expressed by ([!]) were given in [jrjl and |TH] and the special binary Darboux 
transformation was given in Q. 

It is known that the NNV equation system (0) and (0) can be represented as a 
compatibility condition of the linear system 

<$> xy + = 0, (4) 



$t = + + 3v$ a + 3w$ y . (5) 

In general, in order to construct the solutions of a given equation by means of Darboux 
transformation, one may use a fixed solution for a special spectrum parameter of the 
Lax pair. Then one can get a new solution from an old one with Darboux transfor- 
mation. Without spectral parameter in the Lax pair of the NNV equation, we have to 
construct a binary Darboux transformation for the NNV equation. At the same time, 
with the variable separated approach, two arbitrary functions can be entered into the 
solution. 

It is straightforward to see that the NNV system (Q) and @ possesses the following 
trivial solution 

u = 0, v = v Q (x, t), w = w (y, t), (6) 

where Vq(x, t) and Wo(y, t) are arbitrary functions of {x, t} and {y, t} respectively. 

In order to find some new solutions via Darboux transformation and the seed solution 
(ID, the key step is to find the fixed solution of the Lax pair (f|) and (|) with the seed 
©■ 

It is evident that Eq.(f|) have a variable separation solution in the form 



$o = P + q 
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(7) 



where p = p(x, t), q — q(y, t) are two arbitrary functions when u — 0. Substituting Eq. 
(0) and Eq. (§) into Eq. (§) yields 

Pt + Qt= Pxxx + q V yy + 3V P X + 3w q y . (8) 

It is clear that Eq. (§) can be solve by the usual variable separable approach and the 
result reads 

Pt = Pxxx + 3voPx + c{t) , (9) 



qt = q yyy + SwoQy ~ c(t) (10) 

where c(t) is an arbitrary function of t. For given functions Vq, Wq and c(t), it is still 
difficult to solve (^) and (|T0|). However, because of the arbitrariness of the functions vq 
and wo, we can solve the problem in an alternative way. If we consider the functions p 
and q as arbitrary functions, then v and w can be solved out from (||) and (|10|) 

^0 = {?>Px)~ l {Pt ~ Pxxx - C(t)) (11) 



w = (3q y ) (q t -q yyy + c(t)). (12) 

As usual, the binary Darboux transformation for the NNV equation can be con- 
structed by introducing the closed 1-form 

u($, $ ) = ($$ox - $x$o)dx - ($$ 0y - ®o$y)dy + { <I>( <I> , , , - <S> 0yyy ) 

-|-$Q($yyy ^xxx) 2($xx^0x ^x^Oxx ^Oyy^y ^Oy^yy) 

+3v(<S>$ 0x - $ a .$o) + 3w(%% - ®ov$)}dt. (13) 



Then the new wave function is 



$[l] = $ 1 |a; (14) 
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and the equation system Eq. (|j) and @ is covariant with respect to the transformation 
fll4]) and the transformed coefficients u[l], v[l] and w[l] are 

«[1] = 2(ln%) xy =^^- 2 , (15) 

,[1] = v + 2( Mo ) xx = Vo + ^ q -j^, (16) 

w[l] = w + 2(\n%) yy = w + ^L- T ^- 2 , (17) 

p + q (p + g) 2 



where and u>o determined by Eqs.flTTD and flI2|). From (|T5D, we see that the solution 
u[l] obtained by the once Darboux transformation is only a special case of the uni- 
versal expression (HD with 03 = 0. The solution fll5l) -(|17D can also be obtained by the 
truncated Painleve expansion using the similar method as for the AKNS system pl|1 . 

By means of the iteration of the Darboux transformation, we can construct the 
second Darboux transformation from the first Darboux transformation. We also take 
two special solutions of the Lax pair ([|) and ([5]) as the variable separable ones, $1 = 
Pi(x,t) + qi{y,t) and $2 = P2{x,t) + q2(y,t), then the transformed wave function is 

$[2] = J uj{$xM) (18) 

which leads to the formula 

u[2) = 2(ln J u{§ x ^ 2 )) xy , (19) 

v[2] = v + 2{lnju{$ 1 ,$ 2 )) xx , (20) 

w[2) = w + 2(lnJou(<f> 1 ,<f> 2 )) yy (21) 

According to the selection for the arbitrary functions, the potential functions Vq and 
Wq are related to the functions pi, p 2 , qi and q 2 by 

Pit - Pixxx - ci(t) - v p lx = 0, (22) 

qit ~ qiyyy + Ct(t) - W qiy = 0, (23) 
P2t - P2xxx ~ C 2 (t) - V p 2x = 0, (24) 

q2t - q2yyy + C 2 (t) - W q 2y = 0. (25) 
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Now we can consider one of the functions p% = pi(x, t) and p 2 = p 2 (x, t) as an arbitrary 
function of {x, t} and one of the functions qi = q±(y, t) and q 2 = g 2 (j/ ; t) as an arbitrary 
function of {y, t) while the remained functions of v Q , u> , Pi, p 2 , qi and q 2 should be 
determined by (g|), (||), (||) and (H). 

It is not easy to solve the equation system (^), (^3|), (0) and (j25|). However, if we 
choose 

qi = 0, p 2 = 0, , (26) 

then the Eqs. (fl9|)- (|2~T|) can be simplified into 

u[2] = 2\n(c + p 1 q 2 ) xy , (27) 
v[2\ = Vo + 2\n(c + Pl q 2 ) xx , (28) 
w[2] = Wo + 2 ln(c + q\q2) V y, (29) 

where c is an arbitrary integral constant, p\ and q 2 are arbitrary functions of {x, t} 
and {?/, t} while v and ty are fixed by ( fZ2"| ) and (^). It is easy to prove that the 
results (^7|)-(|29|) are equivalent to those obtained via the usual variable separation of 
the multi-linear equation |13|]. Actually, by re-writing c, pi and q 2 as 



CL\d 2 CL\ . . 

c = -a , pi = a 2 + a 3 p, q 2 = q H , (30) 

the negative value of the right hand side of (^) is transformed to the universal quantity 
expressed by ([!]). 

Various types of coherent localized structures for the physical field u of the NNV 



system had been described in [O] and [13 thanks to the arbitrariness of the functions 



p and q. We do not repeated these known localized excitations but write down only one 



new type of localized solutions for the field u expressed by ( |27|) with (|30|) to complement 
the results of [ O . 

In addition to the continuous localized excitations in (l+l)-dimensional nonlinear 
systems, some types of significant weak solutions like the peakons|22"|, |23[| and com- 



pactons have been attract much attention of both mathematicians and physicists. 
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In [[15] and [23|, the possible (2+l)-dimensional localized peakons had been given. The 
(l+l)-dimensional compacton solutions which describes the typical (l+l)-dimensional 
soliton solutions with finite wavelength when the nonlinear dispersion effects were firstly 



given by Rosenau et al.[2^] and may have many interesting properties and possible 
physical applications fl2"5|-f2lj| . For instance, the compacton equations may be used to 
study the motion of ion-acoustic waves and a flow of a two layer liquid[]27|. In [26|, the 
Painleve integrability of two sets of Korteweg-de Vries (KdV) type and modified KdV 
type compacton equations are proved. Because of the entrance of arbitrary functions 
in the (2+l)-dimensional nonlinear physics models, it is easy to find some types of 
multiple compacton solutions by selecting the arbitrary functions appropriately. For 
instance, if we fixed the functions p and q as 

N ( X + dt < X oi - ^7 

p = ^2 { bi sin ( k i( x + c it ~ x 0i)) + k x 0i - afr < x + c4 < x 0i + ^ , (31) 

i=x [ 2bi x + > x 0i + ^ 

and 

m ( y<yoj-i- 

Q = ^2\ dj sm(li(y - y 0j )) + dj y 0j - ^ < y < y 0j + . (32) 
i=i [ 2dj y > y /+ i- 

where bi, ki, x 0i , dj, lj and yoj are all arbitrary constants, then the solution (|7|) with 

(|30D becomes a multi-compacton solution. 

From ( j31~l) and ( |32"|) , one can see that the piecewise functions p and q of the compacton 

solutions are once differentiable 

n ( x + c{t < x 0i - 

Px = ^ { biki COS ( k i( X + ~~ X 0i)) X 0i ~ Wi < X + Cit - X 0i + Wi ' ( 33 ) 

and 



i=l [ X + Cit > X 0i + ^7 



m o y<y 0j -i- 

= { djlj cos ( l i(y - yy)) yoj-i- <y< m + w, > ( 34 ) 
t =x [ v > yoj + w t 

So the multi-compacton solutions expressed by (P (i.e., (|27|) with (|30|)) are still con- 
tinuous any where though there are some isolated non-differentiable lines. 
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In (l+l)-dimensions, a non-differentiable solution, u = Uq, like the the compacton 
(and/or peakon) is called a weak solution of a nonlinear (l+l)-dimensional PDE (par- 
tial differential equation) 

F(u, u x , ut, u xx , ...) = F(u) = (35) 

under the meaning that though the compacton (and/or peakon) solution is non- dif- 
ferentiable at some points x = Xi(t), the distribution, F(u ) = f(S(x — Xi(t)) (where 
5(x — Xi(t)) is a Dirac 5 function and / is a function of the Dirac 5 function and its 
derivatives) is really a zero distribution that means f^°° ftp(x,t) = for arbitrary tp. 

However, in (2+l)-dimensions, the (2+1) -dimensional compactons (and the peakons 
reported in |Tj|) are exact solutions of some (2+l)-dimensional nonlinear equations are 
guaranteed by the variable separation procedure. When we substituting the piecewise 
solutions into the nonlinear PDEs, the Dirac delta function(s) in x (y) directions will 
be vanished by the differential operator in other direction, d y (d x ). In other words, we 
can take (have taken) d y tjj{x,t) = no matter the function ifi is a continuous function 
of {x, t} or a generalized distribution functions (with some Dirac delta functions) of 
{x,t}. 

Fig.l is the evolution plot of a three compacton solution fl27D with (p0|) , (|3l|) , fl32l) 
and 

N = 3, M = 1, ao = 20, ai = a<i = 25a^ = 1, b\ = 63 = — C2 = 03 = —2, 
-b 2 = -a = d x = h = k 2 = fa = h = 1, £01 = x 02 = 2o3 = Vol = 0. (36) 

In [0, we have pointed out that (i) the interaction between two travelling ring shape 



soliton solutions is completely elastic and (ii) the interaction between two travelling 
peakons is not completely elastic, two peakons may completely exchange their shapes. 
From Fig. 1, we see that the interaction between two compactons exhibits a new 
phenomenon. The interaction is nonelastic but two compactons do not completely 
exchange their shapes. 



Fig. 1a 




Fig.2 
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Figure 2. Evolution plot for the quantity u2 = u[2] expressed by ( |39| ) 
at times t = 0. 

Because of the arbitrariness of the functions p and q may have also quite rich struc- 
tures. For instance, (2+l)-dimensional compactons may be only compacted at one 
direction. For conveniece later we call this type of compactons the partial compactons. 
Fig. 2 is a plot of a special partial compacton solution fl2"7|) with (0), (j3l|) 



and 



N — 1, a = 20, a\ = a 2 = 25a 3 = —c\ = k\ = 1, xqi = 0. 



q = 10 tanh(y — t), 



(37) 



(3f 



at t = 0. For the partial compacton shown by Fig. 2, the quantity G possesses the 
following quite simple form (u2 = u[2]) 

\X — t\ > jTY, 

n2 { 100cos(x-t)sech 2 (?/-t) , x _ (39) 



[(90-10sin(x-t)) + (46-4sin(x-t))tanh(y-t)] 2 ' X *' < 2?r ' 
Similar to the first kind of compacton (full compacted) solution as shown in Fig. 1, 

the detailed study of the interaction between partial compactons is also non-completely 

elastic and they will partially exchange their shapes. 

If the arbitrary functions q\ and p 2 are not taken as in (EH) but are solved out from 

(|23|) and (^4[), then we can obtain some further new types of exact solutions which 

can not be obtained by the usual variable separation approach. Actually, when the 

arbitrary functions p\ and q 2 are fixed by ( p2|) and ([25|) at the same time, then the 
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corresponding solution(s) for the functions q\ and p 2 can be found by solving the linear 
equations (^4j) and (|23|). 

In principle, the more kinds of exact solutions can be obtained from further Darboux 
transformations starting from the seed solution {(|18"D, (|T9|), (p0|), (|2~T|)}. However, we 
do not discuss these types of solutions further because of their complexity. 

In summary, the variable separation solutions can be obtained not only by the trun- 
cated Painleve expansion and the generalized multi-linear equations, but also by other 
well known approaches especially by the Darboux transformation. In this short paper, 
the Darboux transformation is successfully used to find the variable separable solutions 
of the (2+l)-dimensional NNV equation. The recursive Darboux transformation may 
yield further new types of variable separable solutions while the new variable separable 
function should satisfy a further constrained condition, say, (|24l) with v being given 

By selecting the arbitrary functions appropriately, one may obtain abundant local- 
ized excitations like the dromions, lumps, ring solitons, breathers, instantons, solitoffs, 
peakons, fractal and chaotic patterns. In addition to these types of localized excita- 
tions, a further type of the localized excitations, compactons, is given in this paper. 
The (2+l)-dimensional compactons discussed here possess the different type of inter- 
action properties as that of the ring solitons and peakons. The interactions among 
compactons are not completely elastic and do not exchange their shapes completely. 
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